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Abstract 

It is shown that the thermodynamic maximum entropy principle predicts negative 
specific heat for a stationary magnetically self-confined current-carrying plasma torus. 
Implications for the magnetic self-confinement of fusion plasma are considered. 
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The goal of the controlled thermonuclear fusion program is to make the energy source 
that powers our Sun available to human society. Deep in the Sun's interior, favorable condi- 
tions for the quasi-stationary nuclear burning of the solar plasma prevail as a result of the 
immense gravitational self-forces that keep this huge accumulation of matter together. Since 
gravitational self-confinement is not operative at the reactor and laboratory scale, alternate 
means of confinement have to be employed to achieve sufficiently high plasma densities and 
temperatures in a reactor. In the perhaps most prominent stationary fusion reactor scheme, 
the tokamak, strong electric ring currents are induced in an electrically neutral plasma to 
achieve axisymmetric magnetic self-confinement in a rotationally invariant toroidal vessel 
T. In a torus with sufficiently large major axis, such a magnetic self-confinement mimics 
gravitational self-confinement on account of the Biot-Savart law, according to which in a 
system of parallel current filaments all filaments attract each other magnetically with the 
same force law as would be the case gravitationally in a system of parallel mass filaments. 
What makes the magnetic forces more attractive (in the double sense of this phrase) than 
gravity for laboratory purposes is their very much bigger coupling constant (~ 10 40 t> 2 /c 2 for 
two electrons moving on parallel trajectories with speed v as measured in the laboratory; 
note that the even stronger (v 2 /c 2 is replaced by 1) electrostatic repulsive forces between the 
same two electrons are very effectively screened in a neutral plasma and are traditionally 
neglected to a good approximation). Of course, the analogy does not extend to all aspects 
of plasma self-confinement. In particular, the solenoidal-vectorial character of stationary 
current densities necessitates the toroidal topology of magnetic self- confinement, whereas 
gravity not only allows but manifestly prefers spherical confinement over toroidal. Unfortu- 
nately, axisymmetric toroidal magnetic self-confinement is not known for its stability either. 
Although major efforts are devoted to the stabilization of the plasma configuration, a vast 
reservoir of instabilities capable of destroying the confinement has dramatically slowed down 
the development of a operating tokamak fusion reactor. 

Matters are not exactly helped by the fact that our theoretical understanding of the 
physics on the various space-time scales that govern magnetic plasma confinement is still 
quite incomplete. In particular, while the solenoidal character of the magnetic induction 
together with the axisymmetry and stationarity of the law of momentum balance tell us that 
the poloidal magnetic flux function \1/ and the toroidal current density j must satisfy some 
local functional relation TZ(^,j,r) = (in which r is the cylindrical distance from the axis 
of symmetry) which turns Ampere's law into some in general nonlinear elliptic partial differ- 
ential equation for \l/ (known in the fusion literature as the Grad-Shafranov equation), the 
actual relation 1Z is not fixed (except for the explicit r-dependence in TV) . Some information 
about 1Z should be contained in the law of energy balance between current drive (through 
an applied toroidal emf and other means), ohmic heating and, ultimately, radiation losses. 
Unfortunately only the so-called classical and neo-classical transport coefficients have been 
computed in some detail jlj whereas the small scale turbulent dissipation mechanisms in a 
tokamak plasma remain a largely challenging open problem. In this situation, theoreticians 
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have been forced to rely on fair judgment and good taste when guessing some additional prin- 
ciple (s) that would effectively complete the characterization of the stationary magnetically 
self-confined plasma torus in a tokamak. 

A sizeable fraction of the literature employs a linear approximation of TZ to get j cx 
rendering the Grad-Shafranov equation linear, and some improved-accuracy modeling uses a 
third-order polynomial approximation of TZ 0. Subsequently the filter of linearized dynam- 
ical stability analysis, based mostly on macroscopic magnetofluid theory and mesoscopic 
kinetic theory, is applied to sort out unstable configurations. While this approach has 
met with a certain limited success, one does not learn what the approximations are ap- 
proximate to. Over the years a number of plasma theorists [[| have argued that an equi- 
librium thermodynamics-inspired maximum entropy principle with a few global dynamical 
constraints^ should give answers close to the truth. In essence the various formulations in 
give for TZ the answer j oc exp(\I/), which leads to a nonlinear Grad-Shafranov equa- 
tion that may have more than one solution depending on the domain T, the boundary 
conditions for and the values of the physical parameters of the problem. In addition, 
this approach provides a global stability criterion within the class of axisymmetric states 
satisfying the same dynamical constraints. Only those solutions that maximize the relevant 
relative entropy functional will be globally stable. 

Since in non-equilibrium statistical mechanics the maximum entropy principle has not 
acquired a status anywhere near as fundamental as in equilibrium statistical mechanics, it is 
mandatory to register some arguments in its favor for the case at hand. Thus, the relation 
j oc exp(^) has been shown to be almost universally singled out also by a truly dissipative 
Fokker-Planck approach to stationary magnetically confined plasma @, [§]. The perhaps most 
compelling reason to give it serious considerations, however, is the successful application of 
the maximum entropy approach to the physically distinct but mathematically quite similar 
problems of stationary planar incompressible flows,0 where the vorticity plays a role closely 
analogous to the current density in the plasma torus, and the strongly magnetized pure 
electron plasma in a circular cylinder,^ where the charge density plays that role. In this 
spirit, we have conducted a thorough investigation of the thermodynamic- type maximum 
entropy approach to the magnetically self-confined stationary plasma torus fl5[. 

1 Since any probability density maximizes the entropy relative to itself, a stationary plasma torus is 
necessarily a maximum entropy configuration under some constraints. What makes the maximum entropy 
proposal non-empty is the insistence on only a few global natural dynamical constraints. 

2 The first qualitative predictions based on statistical mechanics of the Hamiltonian system of N point 
vortices were made in || . The quantitative evaluation began with Q| ; impressive agreement with simulated 
flows is reported in Q . Its mathematical rigorous foundations are by now almost complete, the latest word 
being Q ; see JTo| for a review. More recently a formulation based directly on continuum vorticity has gained 
much ground; see [[n] for a state-of-the-arts report. 

3 In theguiding center approximation the dynamics of this plasma system is identical to that of N point 
vortices [[12j. Statistical mechanics in the corotating frame predicts that at high enough effective energies 
the nonlinear m = 1 diocotron mode has higher entropy than any other configuration with the same energy 
and angular momentum |]l3|] , in accordance with remarkable real experiments fbfl] . 
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A most curious finding of the study is that the gravity-inspired toroidal magnetic 
plasma self-confinement scheme inherits from the stars their gravo-thermal negative specific 
heat!^ This result is a little surprising, for it follows from what we said earlier that the 
plasma torus should actually more closely mimic a cylindrical caricature of a star, and the 
specific heat of a 'cylindrical maximum entropy star' |20j and its plasma physical clone, the 
cylindrical Bennett pinch |2]J , is non- negative! The existence of a maximum entropy plasma 
torus with negative specific heat is therefore a truly nontrivial fact. The purpose of this note 
is to point out some potentially important consequences of our finding for plasma physical 
applications.^ To pave the way for the discussion we first describe the model and our results. 

The model 

Since the whole problem is rotationally invariant, we work with conventional cylindrical 
coordinates r, 8, z. The magnetic induction field decomposes accordingly as B = By + Bp, 
where Bp||e# and Bp _L eg. In an actual tokamak, the toroidal component Bp and a part, Bo, 
of the poloidal component are externally generated harmonic fields that serve the purpose of 
azimuthal and radial stabilization. The total poloidal induction Bp = V^P x V# is the sum 
of B and a component which is generated by the electric plasma current density vector, jeg, 
via the toroidal Ampere's law 

- rV • (r" 2 W) = 4vrc~ 1 j. (1) 

In a similar manner one decomposes the electric field, E = Ep + Ep, where Ep is driving 
the plasma current while the poloidal part is determined by Coulomb's law V • E = inp, 
where p is the electric charge density. Usually the so-called quasi-neutrality approximation 
is invoked, which determines the poloidal electric field in leading order through a singular 
perturbative approach to Coulomb's law. To keep matters as simple as possible, we consider 
an 'electron-positron' plasma, which is totally charge symmetric with regard to the particle 

4 Ref. flUl , pp. 60-63, explains why a homogeneous piece of "everyday matter" must have positive specific 
heat. See p. 62 of the same reference for why those arguments do not rule out negative specific heat in 
an isolated gravitating system. Indeed, the virial and the equipartition theorems imply that in a spherical 
equilibrium system the energy is distributed -2:1 between gravitational and kinetic. A decrease in total energy 
E of a gravitational equilibrium gas ball will increase its thermal energy. Such a system grows hotter while 
losing energy through, say radiation. Negative specific heat in self-gravitational perfect gases is evaluated 
quantitatively already in [p^| and is further discussed in however, none of these configurations with 
negative specific heat is thermodynamically stable though some are metastable. Thermodynamically stable 
self-gravitating configurations with negative specific heat can occur when the Newtonian — 1/r singularity is 
stabilized either as in quantum mechanics or in classical hard balls systems jl9| . 

5 Recently, the existence and importance of negative specific heat was also reported for the diocotron mode 
of the guiding center plasma JD| alias point vortex gas, and for certain vorticity structures in geostrophic 
flows [[llj |2^]. However, very different from the gravo-thermal type negative specific heat that we report here 
to be a characteristic also of the magnetically self-confined plasma torus, the negative specific heat of these 
quasi-particle systems does not couple to the thermal motion of the underlying physical particle systems, 
which is evident from the fact that these quasi-particle systems also exhibit negative temperature jq]. 
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species so that p vanishes exactly. In that case the poloidal electric field vanishes identically, 
too, while the toroidal one is implicitly contained in the electric plasma current / = Nqu/27r. 
Here iV is total number of plasma particles, q the elementary charge, and u the mean absolute 
angular frequency of a species. This settles the electromagnetic part of the model, and we 
turn to the statistical mechanics part. 

While most works in || are formulated at the macroscopic magnetofluid level, they can in 
essence be recovered from the statistical mechanics approach of Kiessling et al. in ||, which 
begins with the Hamiltonian N particle formulation and takes the kinetic limit. At this 
kinetic level, the plasma particles are described by distribution functions on single-particle 
phase space. We seek those distribution functions which maximize the familiar Boltzmann 
entropy functional under the constraints that the two separating integrals of motion, particle 
number and energy, take prescribed values A" and E, and given that the plasma carries a 
prescribed electric current I. Since the current is not a separating integral, one has to resort 
to a ruse and prescribe each species' canonical angular momentum, which in the axisym- 
metric kinetic limit is a separating integral, and subsequently pass from this microcanonical 
angular momenta ensemble to its canonical convex dual, characterized by prescribed u, viz. 
/. The solutions of the corresponding Euler-Lagrange equations for this variational prin- 
ciple are also stationary solutions to the axisymmetric kinetic equations of Vlasov. Over 
velocity space the resulting distribution functions are simply rigidly rotating Maxwellians 
with temperature T = (A^b/?) 1 > 0, tied to the energy constraint, and angular frequen- 
cies ±u; (oc /), microcanonically tied to the angular momentum constraints but canonically 
prescribed. This allows one to explicitly integrate over the velocity space to retain only 
the effective macroscopic entropy principle for the total number density of plasma particles, 
n(x) = n(r,z), which in our charge symmetric plasma is just twice the value of the com- 
mon space-dependent Boltzmann factor of each species' distribution function. In effect this 
renders the entropy a functional of n, 

S(n) = -ka J n(x) In (A| B ri(x)/2) dT + ~Nk B , (2) 

where AdB = hj ^Iwrnk^T is the thermal de Broglie wave length. This entropy functional 
has to be maximized under the constraints that n > is axisymmetric, JndT = N, and 
that the effective energy functionalf] 

W{n) = ~\J J n{x)K(x, x')n(x')dTdT + ^Nk B T (3) 

6 The negative sign in front of the magnetic energy in W is due to the canonical constraint of prescribed 
electric current; cf. the negative sign in front of the centrifugal contribution to the kinetic energy of a rotating 
thermal system in the co-rotating frame, see Landau-Lifshitz (op. cit., pp. 71-73). Incidentally, those very 
centrifugal contributions to W are negligible in our plasma and have been omitted. Moreover, the toroidal 
field Bt does not show since we consider only axisymmetric configurations with toroidal current density. 
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takes a prescribed value, say E. Here, K(x, x') = (2ttI/cN) 2 G(x, x'), and G is the Green's 
function for —V • (r _2 V) in T for boundary conditions detailed below. With the help of 
a variant of Moser's corollary |23j of the Trudinger-Moser inequality it can be shown that 
given J, the entropy functional S(n) takes its finite maximum on the set of non-negative 
axisymmetric densities n(x) = n(r,z) satisfying Jn(x.)dT = N > and W(n) = E, and 
the maximizer is a regular solution of the Euler-Lagrange equation.^ We remark that more 
than one maximizing density function n might exist, and in addition the nonlinear Euler- 
Lagrange equation may have other types of solutions. We call a solution S stable if it is a 
global maximizer of the entropy (for the given constraints), S meta-stable if it is merely 
a local maximizer, and unstable otherwise. Of course, only those S stable solutions which 
exhibit magnetic self-confinement are of interest. 

Results 

Explicitly carrying out the variations and converting the Euler-Lagrange equation for n into 
an equation for using \I/(x) = c _1 J G(x, x.')j('x')/r > dT' and j(x) = qn(x)ujr, we obtain 



Pfirsch's [24] nonlinear Grad-Shafranov equation 



which is to be solved in the torus T for the boundary conditions encoded in G. Solving 
(fD is in general only possible numerically on a computer. However, some explicit analytical 
control is available if one simplifies the actual laboratory geometry somewhat and considers a 
torus T with rectangular cross section {r, 9, z\r,i < r < r Q ;9 fixed; < z < H}. The poloidal 
flux function \I/(x) = ^(r, z) is assumed to satisfy periodic conditions at the z boundary and 
to be constant at r« and r D , so that the radial component of Bp vanishes at the inner and 
outer boundaries of T. In this setting the harmonic poloidal part is simply a homogeneous 
Bo||e 2 , which we choose so that ^(r^z) = ty(r ,z). By gauge freedom we can now even 
set $?(ri,z) = 0. Beside the desired self-confined configurations, these boundary conditions 
allow also unconfined ones, namely Pfirsch's toroidal sheet pinch ||24j| , given by the following 
e 2 -invariant solution of (HI), 



2c ^ cosh(K 2 [2r 2 - r 2 - r 2 ]/2) 
flujq cosh(/€ 2 [r 2 — rf]/2) 



= --Ez in 3i:j %;r ; . (5) 



with k G (0, 00) a parameter and (3(k 2 ) given by 

(3 = Ac 2 q- 2 N- l uj- 2 HK 2 tanh(/t 2 [r 2 - rf]/2) . (6) 

While these solutions do not describe a magnetically self-confined plasma torus, they serve 
as our jumping off point for the numerical computations of the confined configurations. Our 



7 These are quite nontrivial facts. In particular, all this is not true if we relax the condition of axisymmetry. 
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strategy, which was also contemplated by K. Schindler,^ is to look for ringlike bifurcations 
from the toroidal sheet solution (|5]). At an infinite sequence of discrete values E\ > E 2 > 
with E k \ E^ = — 7r 2 / 2 (r 2 — r 2 )/2Hc 2 , other solutions bifurcate off of the sheet pinch 
sequence, breaking its z invariance. The bifurcation points are determined by setting \I/(x) = 
^pf( r ) + €ip(r, z) + 0(e 2 ), with ip(r ita , z) = and ip(r, z + H) = ip(r, z), and expanding 
to first order in e, giving the linearized problem 

- V ■ (r~ 2 W>) + Vi) - V(ip) = 0, (7) 

where (ip) = J ip(r,z)V{r)dT / jV(r)dT, and 

V(r) = -8K 4 sech 2 (/t 2 [2r 2 - r 2 - r 2 ]/2). (8) 

By Fredholm's alternative [p5], the solution of (^) is trivial except for certain discrete values 



of k at which the bifurcations occur. We have proved || that for our T all bifurcations off 
of the sheet pinch are due to modes ipk, k = 1, 2, that satisfy (ipk) = 0. The first mode is 
of the form ip\(r, z) = R(r) cos(27r[2: — zq\/H), with zq arbitrary, and with R(r) satisfying 

-r (r^R')' + (2ir/H) 2 R + r 2 VR = (9) 

for R(ri) = R(r Q ) = 0. With realistic domain dimensions rj = 1, r a = a/2, and H = 2, a 
standard Runge-Kutta solver finds the unique nontrivial solution at k = K\ = 1.62, giving 
E x = 2.72 W„ with W, = l^rj 2 /25c 2 \ Numerical solutions of (|), with n = 1, r Q = a/2, and 
H = 2, were then computed with a well-tested bifurcation code [EfJ], based on a continuation 



method |27j • Our code reproduced the analytical sheet pinch solution and its first bifurcation 
point in excellent agreement with our independently obtained semi-analytical results. We 
then numerically followed the first bifurcating branch that emerges from ipi to nonlinear 
amplitudes, where it develops into a toroidal ring pinch with a double X magnetic structure 
similar to the double X structure in the PDX-PBX Tokamak experiment in Princeton.^ 
Retrospectively, this vindicates our choice of boundary conditions. 

(FIG. 1) 

Our primary interest is in the energy-entropy diagram. Shown in Fig. [| is AS(n) versus 
W(n), where AS{n) = S(n)-Nk B In ^ (Ane/N) 5 (I /Hc) 6 m 3 rf, with e the Euler number, and 
with the density function n running along the computed bifurcation sequences of ring and 
sheet pinch. At sufficiently high effective energies, Pfirsch's sheet pinch is the unique solution 
of (|]) for the stipulated boundary conditions, hence maximizing entropy. Numerically it 
appears to be the case for all W{n) > E\ = 2.72 W,, see Fig. |^. For all W{n) < E\ down to 
W (n) = —0.5 W, where we terminated the computation, the ring pinch has higher entropy 
than the toroidal sheet pinch at the same effective energy. By asymptotic analysis we found 
that also for W{n) \ — oo, and by continuity for W(n) —W», the maximum entropy 

8 Private communication. 

9 For design information, see http://www.pppl.gov/oview/pages/pbxm_design.htm,. 
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configuration consists of a highly concentrated ring pinch which, in rescaled coordinates 
centered at the density maximum, converges to Bennett's cylindrical pinch [21] as W(n) \ 
— oo. On the basis of this evidence we surmise that the ring pinch has maximum entropy 
for all W(n) < Ex, implying its S stability in the class of rotationally invariant plasma with 
effective energy Win) < E\ and current I. We remark that the first bifurcation off of the 
toroidal sheet pinch into the S stable toroidal ring pinch branch is then a symmetry-breaking 
second-order phase transition. 

It remains to determine the specific heat of the configurations, which we recall is neg- 
atively proportional to the second derivative of S with respect to E. Thus we inspect the 
curvature of the graphs of the entropy as function of energy for the various solutions, given 
in Fig. 0. The graph representing the sheet pinch is concave. However, the graph for the ring 
pinch is manifestly convex over the whole computed range of energies —0.5 W, < W(n) < E\. 
We were also able to prove the convexity analytically to second order in perturbation theory 
away from the bifurcation point. This confirms what we have announced earlier: the specific 
heat of the ring pinch is negative! 

(FIG. 2) 

Discussion 

At last, we discuss the potential implications that our finding of gravo-thermal type negative 
specific heat has for the problem of toroidal magnetic self-confinement of plasma. In a 
gravitationally bound plasma, negative specific heat on one hand aids the ignition of nuclear 
burning in a proto-star by heating it up when it loses energy by radiation, but it is also 
responsible for some more spectacular instabilities once the nuclear burning expires, like the 
onset of the red giant structure |T7j] . It would be intriguing enough if the negative specific 
heat of a magnetically self-confined plasma torus should be confirmed to aid the ignition 
of nuclear burning in a tokamak. For this to be so, one would have to be able to hold N 
and / fixed and secure the toroidal invariance (which is what one wants to achieve anyhow), 
while E would have to decrease (the plasma radiation would seem to help in this respect) 
slow enough so that one would essentially evolve along the ring pinch branch in Fig. 2 to the 
left, thereby heating up the plasma while pinching it more strongly. This would not seem 
unwelcome. For now, however, energy leakage by radiation is a serious problem, while at the 
same time the emf current drive leads to yet uncontrolled ohmic heating of the plasma. In 
this case where E is allowed to fluctuate too widely, the negative specific heat will have a very 
unwanted effect on the confinement. This can be illustrated by considering the temperature 
rather than energy E to be controlled by the competition of ohmic heating and radiation 
(still assuming N and / fixed, and toroidal invariance). In that case the canonical ensemble 
determines the stability. But microcanonical and canonical ensembles are not equivalent 



when the microcanonical one exhibits states with negative specific heat [18, 0, Eq, E3], and 



sure enough, none of the computed ring pinches with negative specific heat minimizes the 
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free energy functional 



F{9) = ~ j r-^WfdT+Np-Hn J e^ /c dT 



(10) 



Actually, F is unbounded below for these (3, N and I values, any minimizing sequence 
concentrating on a singular ring current; cf. |2I| for a good discussion of the translation- 



invariant analog. Of course, a real plasma would not get anywhere near such a singular ring 
current configuration, for a highly concentrated plasma ring is known to be susceptable to 
magnetofluid dynamical instabilities that destroy the axisymmetry. 
Conclusion 

To summarize, the S stable ring pinches have negative specific heat of the gravo-thermal 
type and will therefore be stable if and only if N, E, and I are essentially fixed and the 
toroidal invariance is secured, in which case the negative specific heat may aid the ignition 
of thermonuclear burning. The "if" part is good news; the bad news is the "only if" part. 
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Figure 1: Poloidal magnetic lines of force of maximum entropy solutions near the second 
order phase transition at E 1 = 2.72W.. Ring pinch (left): W(n) = 2.34 W„ (3 = 0.29N/W., 
z = 0; Sheet pinch (right): W(n) = 3.00 W„ (3 = 0.30iV/W^.. The toroidal hoop effect is 
neatly visible. 




0.0 1.0 2.0 3.0 4.0 

w/w. 



Figure 2: AS(n) versus W{n) for toroidal sheet pinch (concave branch) and ring pinch 
(convex branch). 
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